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Abstract 

Here /(i?)-cosmology is discussed using a different approach. This 
model explains early-inflation, emergence of cosmic background radi- 
ation at the exit from inflation, cosmic deceleration during radiation- 
dominance followed by deceleration due to curvature- induced matter 
and acceleration in the very late universe due to curvature- induced 
phantom dark energy. This model predicts collapse in the future uni- 
verse. Further, a possible avoidence of collapse as well as revival of 
very early universe is suggested. 

1. Introduction 

Observation of cosmic acceleration in the very late universe [1, 2] 

created lot of sensation in the arena of cosmology and emerged as 

the most fundamental theoretical problem. According to Friedmann 

equations (giving cosmic dynamics), cosmic acceleration is caused due 

to negative pressure p < — p/3 of the dominant fluid (p being the 

energy density). Many field-theoretic, Chaplygin gas [3, and references 
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therein] and /(i?)-dark enenrgy models [4, 5, Q]{R being the scalar 
curvature and f{R) being some function of R) were proposed to explain 
late acceleration and to probe a proper source of the exotic fluid, which 
is invisible but having gravitational effects. This fluid is popularly 
known as dark energy. Here /(i?)-cosmology is adressed. 

In /(i?)-dark energy models, non-linear terms of curvature is identi- 
fied as DE term . But, a more appropriate theory should begin without 
imposing any pre-condition or identifying some terms as DE terms a 
priori. 

This article presents a different approach to /(i?) -cosmology, where 
no curvature term is recognized as gravitational alternative of DE a 
priori[7, 8, 9] contrary to the approach of refs.[4, 5, 6]. 

In refs.[4, 5, 6], gravitational equations are derived from the action 
having Einstein -Hilbert term and non-linear curvature terms. Terms 
in gravitational equations, due to non -linear curvature terms, are rec- 
ognized as DE terms. In what follows, trace of /(i?)— gravitational 
equations are obtained yielding an equation for scalar curvature R. In 
the homogeneous space-time, this equation reduces to the second-order 
equation for the scale factor a{t). First integral of this differential equa- 
tion yields the Priedmann equation giving dynamics of the universe. 
Interestingly, the Priedmann equation, obtained so, contains terms for 
quintessence-like DE term and dark radiation in the early universe as 
well as dark matter and phantom-like DE term in the late universe. 
These terms are induced by curvature without identifying any curva- 
ture term as DE a priori. Here DE terms emerge from Einstein term 
proportional to R as well as non-linear terms of R spontaneously. 
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The gravitational action of the model being addressed here contains 
non- linear terms and R^'^^^^ (with r > being a real number). It is 
interesting to find that dark matter is induced in the present set-up, if 
r = 3. 

The present model yields an interesting cosmological picture from 
the early universe to the future universe. Here, investigations begin 
from the Planck scale (the fundamental scale). It is found that the 
early universe infiates for a short period, driven by curvature-induced 
quintessence dark energy. During this period, elementary particles are 
produced and lot of energy is released as radiation at the end of infla- 
tion. The emitted radiation thermalizes the universe rapidly up to the 
temperature ~ lO^^GeV. As a consequence, produced particles due to 
the decay of curvature-induced quintessence attain the thermal equil- 
librium with radiation. There are two sources of radiation (i) emitted 
radiation at the end of inflation and (ii) curvature-induced radiation 
recognized as dark radiation. The emitted radiation, is identified with 
the cosmic background radiation with very high initial temperature 
~ lO^^GeV. Due to its dominance, the universe decelerates after exit 
from inflation heralding the standard cosmology. 

The produced elementary particles, during inflation, undergo various 
processes of the standard cosmology such as nucleosynthesis, baryosynthesis 
and hydrogen-recombination, which are not discussed here. Thus, like 
radiation, we have two types of pressureless matter (i) baryonic matter, 
which is formed due to nucleosynthesis and baryosynthesis of elemen- 
tary particles produced and (ii) curvature-induced non-baryonic mat- 
ter identified as dark matter. After sufficient expansion of the universe. 
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the matter dominates over radiation causing decelerated expansion as 
~ t^/^ in the late universe. 

It is interesting to note that alongwith curvature-induced phantom 
terms p^[i — {p^)/'^M (Pde being the phantom energy density) in the 
late curvature causes another constant A, which is analogous to negative 
brane-tension in Randall-Sundrum 11 theory of brane-gravity [10]. As 
a remark, it is nice to mention that these type of terms also appear 
in Freidmann equation based on the loop quantum gravity [11]. As 
brane-theory prescriptions are not used here, the curvature-induced A 
appearing in this model is identified as cosmic tension like the Refs. 
[8]. Curvature-induced phantom dominates at the red-shift z = 0.303 
causing a cosmic jerk. As a consequence, a transition from deceleration 
to acceleration takes place in the very late universe. It is found that 
the phenomenon of cosmic acceleration will continue in future too. The 
accelerating phase will end when phantom density will be equal to twice 
of the cosmic tension followed by deceleration due to rc-dominancc of 
matter. 

Here it is shown that phantom density will still grow due to expansion 
of the universe during re-dominance of matter and expansion will stop 
at time tm — 3.45 x IQ-^^^o (^o being the present age of the universe). 
So, the universe will bounce causing contraction in the universe. As 
a result, matter energy density will increase rapidly such that energy 
density and pressure density diverge as well as a = at icoi — 3.62 x 
IO^Hq. This result predicts collapse of the universe at time icoi- 

Further, it is discussed that the result of cosmic collapse is obtained 
using the classical mechanics. It is argued that near tcoi energy density 
and curvature will be very high. It is analogous to the state of the 
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early universe. So, production of quantum particles will take place due 
to rapid change in topology of the space-time in the vicinity of tcoi- It 
is found that back-reaction of produced particles will avoid the cosmic 
collapse and universe will expand exponentially when t > tcoi- 

Natural units {ks — H — c — 1) (where ks, h, c have their usual 
meaning) are used here. GeV is used as a fundamental unit and we 
have IGeV"^ = 6.58 x IQ-^^sec and IGeV = 1.16 x W^^K. 

2. f{R)- gravity and Friedmann equations 

The action is taken as 

S = J d'xV^lj^ + aR' + , (2.1) 

where G = Mp^{Mp = lO^'^GeV is the Planck mass), q; is a dimension- 
less coupling constant, is a constant having dimension (mass)*^^'^''^ (as 
R has mass dimension 2) with r being a positive real number. 
The action (2.1) yields gravitational field equations 

11 1 

-XRn. - 7^9^.uR) + a(2V^V,i? - 2g^,aR - -g^,R^ + 2RR^,) 



+/?(2 + r)(V^V.i?(^+'-) - g,.aR('^^^) + ^(3g,,R^'+^^ 

-/3(2 + r)i?(^+'-)i?^, = 0, (2.2) 

where stands for the covariant derivative. 
Taking trace of (2.2), it is obtained that 
R 



■[6a + 3/?(l + r)(2 + r)i?'']ni?-3/3r(l+r)(2+r)i?('^-^)v''i?V/,i? 

+PrR^^+''^ = (2.3) 

with 
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In (2.3), [6a + 3/3(1 + r)(2 + r)7?''] emerges as a coefficient of OR due 
to presence of terms aR^ and fiR^"^^'^'^ in the action (2.1). If a = 0, 
effect of R^ vanishes and effect of R^'^^'^^ is switched off for (3 — Q. So, 
like [7] an effective scalar curvature R is defined as 

T^*- = [6a + 3/3(1 + r)(2 + r)i?n, (2-5) 

where 7 is a constant having dimension (mass)"^'' being used for di- 
mensional correction. 

Connecting (2.3) and(2.5), it is obtained [9]that 



IGttGjR 



T — 1 



jR"- 



+ nR + (r- l)R-^V''RVi^R 



-iRj'-^ ~ ~ . , . ^ 

-(1 - r)— ^V^i?V^i? + rR-^V^'RV^R 

6a — ^R^ 

Experimental evidences [12] support spatially homogeneous and flat 
model of the universe 

dS'' ^ - a''{t)[dx'' + dy" + dz"] (2.7) 

with a{t) being the scale factor. 

For a(t), being the power-law function of cosmic time, R ~ a~". For 
example, R ~ for matter-dominated model. So, there is no harm 
in taking 

^=4' (2-8) 

a"- 

where n > is a real number and Ais constant with mass dimension 
2. 

Using (2.6), (2.6) is obtained as 



-+ 
a 



, _ , nil — r^-yA^a 
2—n—n[r—l)-\ — nr 



6a - 7A^a-"^ WaJ 1671^7^^ L 7A' 



-1 
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xor . X. xni^i/-[6a-7A'"a-"n2+V^ (2.9) 
n(7A^a-"^)2[3r(l + r)(2 + r)]i+V'^^ / J V y 

Approximating (2.9) for small a{t) in the early universe and inte- 
grating, we obtain 

a\2_ B 2/3-1/'- 



(2.10) 



a J a(2+2M) n(7^'-)2[3r(l + r)(2 + r)]i+V^a(2+2M) 

with S being the integration constant and 

M ^ 2 - n - nr. 



(2.11) 



Approximation of (2.10)for large a(t) in the late universe leads to 



a 

- + 
a 



2 - 2nr] = Da""' - Ea^""' , 



(2.12a) 



where 



D 



(2.126) 



and 



/^_^Vr^^ [3r(l + r)(2 + r)]-i-V- /6axi 



.7^^/ Ll67rGn 
(2.12a) is integrated to 
2 B 2D 



(-) =- 



{2+2N) {2 + 2N + nr) 



E{2 + 2N + nr) 
D{2 + 2N + 2nr) 



with 



N ^2- 2nr. 



(2.13) 



(2.14) 



Further, it is found that if M = 1, the first term on r.h.s. (right hand 
side) of (2.10) gives radiation. Moreover, the first term of r.h.s. of 
(2.13) has the form of matter density if A?" = 1/2 . So, setting M — 1 
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in (2.10) and = 1/2 in (2.14) to get a viable cosmology, it is obtained 



that 



and 



n — 



r = 3. 



(2.15a) 



(2.156) 



3. Power-law inflation, origin of matter and reheating after inflation 
3(a). Power-law inflation The Priedmann equation (2.10) for the 



early universe is approximated to 



^4+ 3 ^de 



using (2.11) and M = 1 for 

/ \ 4/3 

a < ( 7y4^/6aj = a. 

Using (3.2), the energy density in (3.1)is obtained as 

16/3-V3 



Pde — 



SttG 



,3/2 



-3/4 



3/4 



7/3 



(3.1) 



(3.2) 



(3.3) 



.ll(7A3)-i/3[i80]V3 
is caused by linear as well as non-linear terms of curvature. Due to its 
origin from curvature, p^" is identified as dark energy density. 

Interestingly, a radiation density term B/a'^ emerges in (3.1a) spon- 
taneously from gravity. This type of term emerged first in brane-gravity 
inspired Priedmann equation. So, analogous to brane-gravity, here also 
B/a"^ is called dark radiation. 

Here investigations start at the Planck scale, where DE density is 
obtained around W^GeV^. Using p^l = lO^^GeV^ at a = ap, (3.3) is 
re-written as 



qu 



ap/ 



a-3/4_«-3/4^7/3 



-3/4 



dp — Clc 



3/4 



(3.4) 
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(3.4) and the conservation equation 



Pde + 3-(pde+Pde) =0 

a 



yield 



Pde 



7/3 



(3.5) 



(3.6) 



(3.6) shows violation of SEC for ap < a{t) < Og. It implies dark energy 
given by (3.2) mimics quintessence. 



As Op is expected to be extremely small, so dominates over the 
radiation in (3.1). Moreover, (3.6) shows that p^" — at a — Qc- So, 

for ap < a{t) < cosmic dynamics is given by 



a/ 



)■ 

Op/ 



-3/4 



— a 



-3/4 



which integrates to 



giving power-law inflation during the time period 

tc-tp^ 7.77 X 10%, 

if 



/ a \ -1/4 
Vop/ 

(3.7) 



(3.8) 



10 



28 



ap 



(3.9) 



(3.10) 



being required for sufficient inflation. 

3(b). Origin of matter and reheating after inflation 

Ref. [9] demonstrates in detail that the curvature-induced quintes- 
sence DE with density (3.4) can be reahsed through the quintessence 
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scalar (f)o{t) with potential 



V{<f>o) 



-F[a- 



3/4 _ -3/4.4/3 



3/4 _ ^-3/4n _ Z_ 
" ' 30J 



where 



and 



4 ^ 



-3/4 



-3/4 



-7/3 



(3.11a) 
(3.116) 
(3.11c) 



The curvature-inspired 0o(O is the background field deriving infla- 
tion and (j){t^ x) playing the role of inflaton can be realized as (/)(t, x) = 
0o(^) + ^(pit^^) with d(f)(t,:K) being the quantum fluctuation. Here, 
perturbations in the metric components are ignored for simplicity. The 
curvature- induced inflaton (p{t,x) decays into bosons and fermions due 
to fast topological changes during inflation [9, for details] . 

Moreover, fluctuations S(po{t) around (p — (j){tc) given by the equation 



a 1 



(X 

(3.12) yields the solution 

5(l)o{t) ~ ^j2/Tdyq-^'^cos{hr} - 7r/2 + 237r/8) 

with 



77 = 



8V37r 



(3.12a) 



(3.126) 



(3.12c) 



(3.12a) shows the release of energy as radiation at the end of inflation 

due to fluctuations with decaying amplitude. Density of the released 
energy at the exit of the universe from inflation is obtained as 

y(0) - V{(t)c) ^ lO^^GeV^ (3.13) 
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The emitted radiation, having energy density (3.13), reheats the uni- 
verse rapidly up to the temperature 



As a result, created elementary particle are highly relativistic and 
have thermal equillibrium with the emitted radiation. 

4. Deceleration and acceleration in late and future universe 

4(a) .Deceleration driven by radiation The early universe from 

inflation at t = tc, when a = Oc and = as it is given by 
(3.4). As discussed above, the emitted radiation at this epoch re- 
heats the universe upto very high temperature given by (3.14) and 
radiation-dominated era of the standard model of cosmology is recov- 
ered. So, the emitted radiation is identified as the cosmic background 
radiation{CMB) . 

Thus, we have two sources of radiation (i) CMB and (ii) curvature- 
induced dark radiation given by (3.1a). As temperature of CMB, ob- 
tained here, is very high, dark radiation too will have thermal equillib- 
rium with CMB. So, energy density of created particles, dark radiation 
and CMB together will have energy density 



So, at the end of infiation (a = ac),Friedmann equation (3.1a) re- 
duces to 



T; = 4.8 X lO^^GeV. 



(3.14) 




(4.1) 




(4.2) 



and it is integrated to 



a 



(t) = ae[l + AMpy/n/ldatit - ^J^^. 



(4.3) 
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4(b). Matter-dominance and deceleration Like radiation, we have 
two types of matter also (i) dark matter given as g^Q^a on setting 

= 1/2 in (2.13), which is non-baryonic due to its origin from gravity 
and (ii) baryonic matter formed by elementary particles ( produced 
during inflation) through various processes of standard cosmology such 
as nucleosynthesis, baryosynthesis and recombination of hydrogen not 
being discussed here. 

It is interesting to see that if 



and 



(2.13) looks like 

(°)' = ^ + 0.73g„V/«{l-°y| (4.6) 
\a/ a-^ I 2A J 

using current values of matter density p^^^ — Q-Hpq and dark energy 

density pg^ = O.TSpJj' provided by WMAP [13]. Here p"^ = |§ with 

Ho = lOOhkm/Mpcsecond = 2.32 x IQ-^^hGeV = [0.96to]"^ with to = 

13.7Gyr = 6.6 x lO^^GeV"^ and h = 0.68. 

The conservation equation (3.5) for p^^ yields 



w£ = -l (4.7) 



It means that the curvature-induced energy density p^ mimics phan- 
tom . Thus, in the late universe, a phantom model is obtained from 
curvature without using any other source of exotic matter alongwith 
cosmic tension A(mentioned in the introduction). 

Prom (4.6), we find that 0.27 H^/a^ > O-TSH^a^^"^ for a < a* and 
0.21 Hi /a^ < 0.73Hla^^^ for a > a*. It means that a transition takes 
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place at 



/ 23 \ 4/15 
a = a, = (- j = 0.767. 



(4.7) 



It shows that phantom terms in (4.6) dominates over matter term at 
red-shift 



= 1 = 0.303, 



(4.8) 



which is very closed to lower limit of given by 16 Type supernova 
observations [2]. 

For a < 0.767, [0.73p^'a^/^]'^ « 0.73p^'a^/^. So, (4.6) is approxi- 
mated as 

2 0.27H^ 



a) 



which integrates to 



a{t) = a,[l + -V^H^a-/'\t - ta)f/\ 



(4.9) 



(4.10) 



Here, td —— 386kyr = 2.8 x 10 ^to is the decoupling time ( decoupling 
of matter from radiation) and the scale factor aa StX, t — td is given 
by l/ttd ^ 1 + Zd ^ 1090 (WMAP results). (4.10) shows deceleration 
during matter-dominance. 

4(c) .Late acceleration during phantom dominance For a > 0.735, 
(4.6) is approximated as 



-3/4 _ 



0.73p 



2A 



}] 



(4.11) 



with Hq given above 
(4.11) integrates to 

a(t) -- 



2A iV ■ 2A 

3 ; ^ 2n -4/3 

--HoV0?f3{t-Qj J 



(4.12) 
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as a^^^^ = 1.22. (4.12)shows acceleration and it is singularity-free. 

5. Re-dominance of matter, cosmic collapse and its avoidance 
5(a). Deceleration in future universe and cosmic collapse 

(4.6)s1lows that phantom-dominance will end at a{te) = Oe = [0.73p°j,/2A]~^/^ 
ie. = 2Aand matter will re-dominate. So, (4.6) will again reduce to 
(4.9) and the future universe will decelerate as 

a{t) = ae[l + ^]Vo:27Hoaf/''{t - te)f^ (5.1) 

It is notable that phantom density increases as universe expands, so 
> 2 A when a > Ue- As a consequence, p^{l — P^/'^^} < using 
(4.4)in (4.6). As growth of will continue, at a certain value > 
of a{t) 

i^ = 0.73«^{^^-l}. (6.2) 

(5.2) shows that d = at a = Um- So, universe will bounce back and 

contract. As a consequece, dominance of matter will increase and FE 
will have the form of (4.9) yielding 

= (^) ~ -VMlHoa-'/^. (5.3a) 

(5.3a) integrates to 



a{t) = am[l - -^/^Hoa;^'\t - tm)f'^ (5.36) 

showing decelerated contraction as d < 0. 
(5.3b) yields a{t) = at 

t = tool = tm + — i=i/o"'«m^' = 3.62 X WHo. (5.4) 
3v 0.27 

So, at t — tcoi, dominating energy density term 

pimat) ^ o.27p^Va^ (5.5) 
in (4.9), will be infinite. These results show cosmic collapse a,t t — icoi- 
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5(b). Avoidance of cosmic collapse Near the collapse time t = tcoi, 
energy density will be very high . So, like curvature-induced quintes- 
sence scalar in sections, the energy density (5.5) can be realised through 
another scalar x) . Due to decay of x in rapidly changing space- 
time, particles will be created near the collapse time. Energy density 
of created particles is obtained as [9, for details] 

Pcreated = ^V^TT^a^/^ (^) e^'^^.rf(37rV7/2) (t^O-''/'^-'- 

(5.6) 

It is natural to think that created particles will effect cosmic dynam- 
ics. As a consequence, FE (4.9) is modified as 

Stt 

~ 0.27H^a-^ + y Mp ^created (5.7) 

The solution of (5.7) can be taken as 

a = acoiea;p[|D(ieoi - ^)}| + 7l^(^coi - t)\% (5.8) 

where Ocoi — a{t — tcoi) , Disa constant of mass dimension and constant 
7 is dimensionless. 
If (5.8) satisfies (5.7), 

^M^D^ = ^-^MlH^.a-:, + Aa-^, (5.9a) 
^Mj>D'^ = -3^M|ifo'«c~oi - 5Xa-^, (5.96) 



where 

A = ^v^7r^a:/^(^)e^'^^.rf(37r^/2)(r,0-^'/^- 

(5.9a) and (5.9b) yield 7 = — 1|. At the largest energy mass scale 
i.e. the Planck mass, energy density is obtained as Mp/Sn'^. Using it 
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0-81 tt2 -.s , r r , 25-1 
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in (5.9a,b,c),acoi = 2.25 x 10^"^^ and D — ^= are obtained after some 
manipulations. 



6. Salient features and concluding remEirks 



Here, a cosmological picture is obtained from the gravitational ac- 
tion containing the linear Einstein term as well as non-linear terms 
and using an approach different from the work [4, 5, 6]. This 
approach has an advantage to have power to explain (i) power- law in- 
flation in the early universe and graceful exit from this phase, (ii) cre- 
ation of SM particles, (iii) recovery of the standard cosmology with 
the cosmic background radiation with extremely high initial temper- 
ature ~ lO^^GeV, (iv) deceleration of the universe driven by emitted 
radiation during the inflationary phase and particles in thermal equilib- 
rium with radiation, (v) deceleration driven by curvature induced dark 
matter and baryonic matter caused by various processes like nucleosyn- 
thesis, baryosynthesis, hydrogen re-combination of elementary particles 
created during inflation, (vi) dominance of curvature-induced phantom 
at red-shift z — 0.303 (which is consistent with observational results), 

(vii) transient acceleration driven by phantom in the very late universe, 

(viii) re-dominance of matter, contraction of the universe, collapse of 
the universe at time tco\ = 3.62 x lO^^to , (ix) avoidance of collapse due 
to creation of particles near the time tcoi as well as its back-reaction 
and (x) rebirth of the universe after tco\- 

Thus, this model predicts the possible revival of the state of the early 
universe at time icoi = 3.62 x lO^^^o ■ 
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